Study of size effects in the structural transformations in bcc Zr films 
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The effect of the film thickness and temperature on structural phase transformations in infinite 
films of bcc zirconium with (001) crystallographic orientation is studied by the molecular dynamics 
method with a many-body potential of interatomic interaction obtained within the embedded atom 
model. It is shown that the mechanism, sequence and final structures of phase transitions in bcc 
Zr films essentially depend on the film thickness. The films with (001) surface up to 6.1 nm thick 
experience an orientational transition into the bcc (110) phase through an intermediate metastable 
fee phase, and then, on cooling, a diffuse bcc (110) — > hep transition is observed. In films 6.1 to 8.2 
nm thick there forms, on cooling, a twin fee phase as a result of shear deformation, so that the film 
surface acquires stepped relief. With further increase of the bcc (001) film thickness there occurs 
martensitic transformation with the formation of a twin hep structure, and the film has a wavy 
surface. The bulk and elastic moduli are calculated for the bcc and fee lattices. It is shown that for 
the fee phase the lattice stability conditions are satisfied in both bulk and film systems. 



The possibility of the formation of anomalous phases 
in thin films of transition d-metals and the conditions for 
their stabilization have long been discussed in the litera- 
ture [lH3j • Such anomalous phases include polymorphous 
structures not observed in bulk samples, but stabilized 
in thin films. It is well known, for example, that on the 
equilibrium P-T phase diagram for the transition metals 
Ti, Zr, Hf in macroscopic state the fee phase is lacking. 
At low temperature and pressure all of these three met- 
als have a hep structure (a phase) which with increasing 
temperature turns into a j3 phase with bcc lattice. The 
high-pressure bcc phase is also realized at room temper- 
ature. At low temperatures an increase in pressure leads 
to a sequence of structural transformations: hep — > cj — ► 
bcc. In zirconium, for instance, the uj phase exists in a 
pressure range from 2.2 to 30-35 GPa. At the same time, 
it has been experimentally found that in Ti thin films 
grown on NaCl, Al, Ni, SiC substrates a fee structure is 
realized [8 |4lll| . Fee hafnium was detected in Hf/Fe mul- 
tilayers [l2j. In Ref. Q it is stated that fee Zr films thick 
500A were obtained by deposition on glass substrates at 
250 — 400°C, with increasing temperature above 450°C 
a mixture of fee and hep structures was observed. In 
another paper [5] fee zirconium films of 50 monolayers 
were reported to be obtained in the temperature range 
1300K - 1800K. In Ref.[13j were obtained nanoparticles 
of polycrystalline fee zirconium with fee grains of 5-10 nm 
in diameter and a lattice parameter ranging from 0.4686 
nm to 0.4714 nm which was shown to increase with de- 
creasing grain size. However, in another experimental 
work [TJ the anomalous fee phases of Ti, Zr, and Hf are 
considered as interstitial impurity phases of residual oxy- 
gen atoms. The authors cast some doubt on the possi- 
bility of the polymorphous transition realization through 
the mechanism suggested in Ref.p}, adducing thermody- 
namic estimates of the critical film thickness for a phase 
transition to occur in a free polycrystalline film. They 
believe the critical thickness of an hep - fee transition to 



be 6.1 nm. Based on an analysis of their own experimen- 
tal data and those available in the literature, the authors 
of Ref. 0| conclude that "by now (2009) there has been 
no one direct experiment which confirms the formation 
of polymorphous modifications not known in bulk sam- 
ples", and "in ultrahigh vacuum the formation of only 
those modifications occurs which are proper to substance 
in a bulk state". Thus, various groups of experimenters 
disagree on the nature of anomalous phases arising in 
thin films. Some state that the anomalous size-induced 
phases observed are of polymorphous nature, others con- 
sider them as being of impurity character. 

Theoretical calculations of the total energy of the Zr 
ground state within the density functional theory show 
the fee structure to be energetically preferable to the bcc 
one. In Ref. on the basis of first-principles calcu- 
lations of the Gibbs energy for free nanoparticles of ti- 
tanium and zirconium, the critical particle size was de- 
termined as a function of the temperature at which the 
fee lattice becomes stable. At T = 300K this size for 
Ti and Zr is close to 6 nm, which agrees well with the 
experimental data [liij ). 

A possibility for the existence of a metastable fee phase 
in zirconium has also been discussed in some works based 
on the use of the molecular dynamics method. In Refs. 
Til 1||, for example, the fee structure was found to form 



at twin boundaries in the hep phase of zirconium, and 
in Refs. [2(| 21 1 the formation of a fee ZR phase oc- 



curred near the crack tip on the boundaries of twins with 
hep structure. A molecular-dynamics study of marten- 
site structures in Zr nanotubes [22| also revealed regions 
of fee phase connecting the nearest hep domains. Up to 
now, there has been no systematic study of the phase 
size-effect, the conditions of structural stabilization, new 
phase formation and morphology, mechanisms of phase 
transformation in nanomaterials of the Zr group. The 
aim of this work is to investigate the temperature stabi- 
lization region of different crystalline structures in zirco- 
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nium nanofilms depending on their thickness at a con- 
stant pressure (P=0). 



CALCULATION METHOD 

The structure stability and physical characteristics of 
zirconium films were studied by the molecular dynam- 
ics (MD) method using the standard XMD package [24|. 
The many-body potential (2|| constructed within the 
"embedded atom" model (EAM) (2(| was chosen to de- 
scribe interatomic interaction in zirconium. In Ref.(25| 
it is shown that this potential allows one to obtain, to a 
high degree of accuracy, the bcc and hep lattice parame- 
ters of zirconium, cohesive energy, elastic constants, melt- 
ing temperature, and other physical characteristics. The 
bulk phonon dispersion curves calculated with this po- 
tential along the symmetrical directions of the Brillouin 
zone of bcc zirconium are shown in FigJTJ The spectrum 
was calculated with a Fourier transform of the function 
of atomic displacement time evolution averaged over an 
interval of 40 ps. The polarization vectors ek with wave 
vector k were calculated in the harmonic approximation 
as the eigenvalues of dynamical matrix for a bcc lattice 
with previously determined parameters corresponding to 
specified temperatures[27| . A comparison with the exper- 
imental data 28j shows that the potential chosen allows 
one to reproduce the experimentally observed features of 
the Zr phonon spectrum and their variation with temper- 
ature, including the softening of the transverse N-phonon 
with decreasing temperature, and so it may be used in 
simulating the thermal stability of bcc zirconium films. 
In this work, the molecular dynamics calculations were 
performed in two ways: (1) the formation of a crystal- 
lite which then is hold, in the state of free evolution, at 
a specified constant temperature and constant pressure 
(P=0) for a time no shorter 

than t — 100 ps, and (2) the formation of a crystallite 
at a high temperature, holding for a time t = 100 ps, 
and then a slow cooling at a constant zero pressure with 
a temperature step AT and a holding time no lesser than 
t — 50 ps at each new temperature. The time step was 
constant and equal to At = If s. In what follows the first 
way will be referred to as "isothermal holding", and the 
second one as "slow cooling". Free boundary conditions 
were specified along the z axis. The film thickness on the 
z axis varied from N z = 5 to N z = 30 bcc unit cells. In 
the x and y directions cyclic conditions were used with 
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N x>y (N x = N y ) ranging from 7 to 36 u.c. Thus we 
simulated thin films infinite in x, y, with free surfaces on 
the z axis. 

In all variants, the calculation started with the forma- 
tion of a crystallite with perfect bcc structure and (001) 
crystallographic orientation. The lattice parameter was 
chosen equal to its value in bulk bcc zirconium at the 
corresponding temperature [27]. The state of the sys- 
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FIG. 1: Phonon dispersion curves along the main direc- 
tions of the Brillouin zone of bcc Zr obtained at = 300 K 
and T = 900 K. Black circles indicate the experimental 
frequencies [28]. 



tern was followed up using the atomic radial distribution 
function (ARDF), the variation of the total energy, and 
visualization of both the crystallite as a whole and any 
separate crystallite plane. The belonging of each atom in 
the system to one of the following structures: bcc, hep, 
fee, and ui was deduced from its nearest environment. 



SIZE DEPENDENCE OF STRUCTURAL 
TRANSFORMATIONS IN BCC (001) 
ZIRCONIUM FILMS 

Structural transformations under isothermal holding 

The calculations by the isothermal holding technique 
have shown that for the films with (001) surface there is 
a critical thickness N z — 17 u.c. (about 6.1nm) at which 
the sequence of structural transformations and their 
mechanism substantially change. In what follows the 
films with thicknesses N z < 17 u.c. and N z 17 u.c. will 
be referred to as "thin" and "thick", respectively. The 
MD simulation has revealed that bcc (001) thin films are 
unstable over the whole temperature range considered; 
they undergo structural transformations the sequence of 
which is temperature dependent. Under isothermal hold- 
ing in the range 500.fr < T ^ 1300.fi the sequence of 
structural transitions is the following: /3 — > 7 — > (3' — > a. 
Here /3 is the original film with bcc structure and (001) 
free surface; 7 is an intermediate metastable fee phase; /3' 
is a twin bcc structure with (110) free surface, and a is 
a hexagonal phase. Besides, the modification of the final 
phase depends on whether N x>y (the number of cells in 
the basic crystallite along the x, y axes) is even or odd: at 
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FIG. 2: The T — N z phase diagram obtained by "isothermal 
holding" for thick bcc (110) Zr films. is the bcc(llO), and 
P is the original bcc(001) phase, jtw/a is the twin fee phase 
with thin hep interlayers, atw/j is the twin hep phase with 
small fee inclusions. 



odd values of N x y < 20 there forms an orthorhombic a" 
phase, while at any even N x>y an a phase with ideal HCP 
lattice arises. In our opinion the a" phase formation is 
due to the lack of vibrations with wave vectors located 
near the point N (k = 1/2 [110]) of the Brillouin zone of 
BCC lattice, which are necessary for the — > a tran- 
sition to be completed. This point is discussed in more 
detail in our previous paper [29j . At temperatures of thin 
film formation T 500 K, the system energy is not suf- 
ficient for the 7 — > f3' transition which proceeds with an 
increase in energy, therefore the 7 phase remains stable. 
At temperatures above 1300 K the twin bcc structure 
with (110) surface is stabilized. 

In thick films the dependence of structural transforma- 
tions on the crystallite size along the x,y axes is not so 
significant, and the calculated results may be presented 
on the Temperature - Film thickness (T — N z ) phase dia- 
gram shown in Fig[21 The diagram displays the areas of 
final stable phases observed in thick films. Note that the 
orthorhombic phase did not arise in thick films because of 
the hep phase being formed here not through the phonon 
mechanism (cooperative atomic displacement), but as a 
result of shear deformation. In addition, in thick films of 
zirconium the sequence of phase transformations, as well 
as the temperature ranges of the final phase stability es- 
sentially differ from the results obtained for thin films. In 
thick films reorientation of the original bcc (001) lattice 
into the bcc (110) one takes place only at temperatures 
T 1400K. However there is a sufficiently wide tem- 
perature range, increasing with the film thickness from 
[900 K - 1300 K] to [400 K - 1800 K], in which the initial 
bcc (001) structure remains stable. Simulating bcc (001) 
thick films by the isothermal holding technique shows 
that a structural transition from the bcc to a twin fee 




FIG. 3: View of the basic crystallite with N z = Ylu.c, N XtV = 
24u.c. from the [100] direction after isothermal holding at 
T=400 K. 



phase with a change in the film surface relief occurs in 
the range of 400 to 700 K. The bcc structure transforms 
into a twin plate structure where the fee plates are sep- 
arated by thin hep interlayers. The emergence of twins 
on the free surface of the film results in the formation of 
a stepped relief, as shown in Fig|3J The figure demon- 
strates a side view of the (001) surface of a film with the 
basic crystallite dimensions N z = 17 u.c.,N XtV = 24 u.c, 
at a temperature T — A00K. One can clearly see steps 
with terraces comprising a few atoms which arise on the 
surface as a result of shear deformation of fee layers. In 
thick films the number of such twin plates increases with 
N z , and consequently increases the portion of the hep 
phase. 

An increase in the density of fee twin plates is also 
observed with decreasing temperature, i.e. the lower is 
the temperature of isothermal holding with constant ba- 
sic crystallite dimensions N z and N x>y , the larger is the 
portion of hep phase observed in the final twin struc- 
ture of a thick film. Moreover, at T — 200K this kind 
of twinning may result in the formation of a homoge- 
neous hep phase which can be considered as a limiting 
case where each fee plate involves only two atomic lay- 
ers, all plates being shifted relative to each other in the 
[112] direction. Conversely, the higher is the tempera- 
ture, the smaller is the number of fee shearing twins in 
the crystallite, and the smaller is the hep phase portion. 
At certain temperatures, ranging from 500K to 600-ftT, in 
thick films with a thickness 17 u.c. sC N z sC 19 u.c. we ob- 
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served the formation of a homogeneous foe phase which 
remained stable throughout the whole time of simulation. 
With increasing film thickness the temperature range of 
the transition to the fee phase decreases, and also the 
portion of fee phase in the final structure gradually de- 
creases. In films with a thickness 17 u.c. ^ N z ^ 23 u.c. 
the portion of the fee phase in the final structure is larger 
than that of the hep phase, while in the thickness range 
24 u.c. ^ N z ^ 31 u.c. the portion of the fee phase is sig- 
nificantly less than the hep fraction. From N z — 32 u.c. 
and on, the fee phase is not observed at any tempera- 
tures, and only an a phase forms as the result of transfor- 
mation. Thus, the isothermal holding simulation shows 
that in thick films of bec zirconium the phase with fee 
lattice arises only in a narrow range of thicknesses and 
temperatures. Recall that in thin films the fee phase is 
stabilized at temperatures below 500 K for any value of 

So we came to the conclusion that using the isothermal 
holding technique which simulates "rapid quenching" of 
the material, one can obtain a homogeneous fee film no 
thicker than a critical value of 17 u.c, that is, about 
6.1nm. This is precisely the critical film thickness of the 
fee — >• hep transition reported in Ref.[7], which was ob- 
tained using the phenomenological theory proposed in 
Ref. [lj. The formation of a twin (fec/hep) film was ob- 
served in our simulation up to a thickness of 30 — 32 u.c. 
corresponding to about 11 — 12 nm. 
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FIG. 4: Dependence of the bec— >hcp transition temperature 
on the thickness of bec (001) Zr film, obtained by "slow cool- 
ing". N x . y = 24u.c 1 is the one-stage hep phase, 2 is the 
"diffuse" transition to the hep phase, 3 is the fec/hep mix- 
ture formed by the shear deformation mechanism, 4 is the 
twin hep phase formed by martensitic mechanism, with the 
twin orientation in accordance with long-wave vibrations. 



Structural transformations on slow cooling 



We have also studied the phase transformations in bec 
(001) zirconium films by the "slow cooling" technique. 
Thin films were cooled from a temperature T — 1500A, 
and the thick ones from T = 1000A", in both cases with 
a step AT — 10K and a holding time t = 50ps at each 
new temperature. The calculation results are shown in 
FigH] Here the film thickness N z in unit cells is plotted 
on the horizontal axis, the crystallite dimensions along 
the x, y axes being the same in all computation vari- 
ants: N x>y = 24 u.c; the temperature is plotted on 
the vertical axis. The temperatures of different struc- 
tural transformations are marked by different symbols. 
The films with N z = 5,6u.c. after reorientation to the 
bec (110) phase, completely transform at high tempera- 
ture into the hep phase. The temperature of this tran- 
sition is shown by circles. In thin films with a thickness 
N z — 7 — 16m. c. the transition from the bec to hep phase 
does not occur at one temperature value. Instead we 
observe the so-called "diffuse" transition which is char- 
acteristic of martensitic transformations when the phase 
state changes not abruptly but gradually over a certain 
temperature range [23]. In our calculations we observed 
phase transitions occurring in two stages. The initial 
and final temperatures of such extended in T transfor- 
mations are plotted by squares. For instance, in a film 
with N z — 7u.c. the transition began at T — 1340K, at 
that moment a part of the basic crystallite transformed 
into the hep phase, which was accompanied by a small 
jump in the total energy curve. Only when the temper- 
ature was lowered to T = 660A, the remaining part of 
the crystallite changed to the hep structure, again with a 
jump in the total energy curve. The sum of the energies 
of these two jumps was equal to the value of the energy 
barrier between the bec and hep phases obtained in the 
case of the one-stage phase transition. The diffuse tran- 
sition is displayed in Fig [51 where a film of dimensions 
N z — 7 u.c, N x _ y — 24 u.c. is pictured from the [001] di- 
rection at different instants of cooling. The correspond- 
ing temperatures are indicated under each picture. It is 
of interest that the crystallite part which has transformed 
into the hep phase at the initial temperature, chaotically 
moves through the crystallite along the [100] direction, 
slightly increasing in size, and when a critical volume is 
attained, at the final temperature all the remaining part 
of the bec (110) phase abruptly changes to the hep struc- 
ture. Note that if the crystallite is held at an intermedi- 
ate temperature, say, at T = 1000A for 1000 ps (1 mil- 
lion steps), the transition is not completed and the film 
part transformed to the hep phase continuously travels 
through the crystallite without growing in size. A similar 
effect, i.e., the existence in a wide temperature range of 
an intermediate configuration of the uncompleted phase 
transition, was observed in studying the influence of the 
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size of free nanocrystals on the j3 — > a transformation 
in zirconium f2l\). As seen in FigU with increasing film 
thickness the temperature interval between the transition 
onset and completion decreases . For instance, in a film 
with N z — 15 u.c, the transition starts at T = 920A 
and is completed at T = 830A , while for a thickness 
N z — 16 u.c. the transformation into the hep phase oc- 
curs immediately in the whole crystallite at a tempera- 
ture T = 920AT. 
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FIG. 5: The diffuse transition. 

A considerable difference in transition temperature is 
observed in the figure between the thin and the thick 
films, which is due to the fact that in thin films the 
transition proceeds from the reoriented twin bee (110) 
phase, while in thick films with N z = 17 u.c. and more 
it starts from the bec (001) phase. Moreover, the films 
with a thickness N z = 17 — 24 u.c. transform not in the 
hep phase, as might be expected, but in the twin fee 
phase. This phase is marked by crosses in Fig0J Note 
that on cooling these films down to T = 100A" with the 
same procedure their structure does not change, i.e. the 
twin fee phase remains stable. With an increase in the 
film thickness the temperature of the bec— sdec transfor- 
mation slightly decreases. In the arising twin fee phase, 
thin interlayers of hep phase are observed along the twin 



boundaries. As in calculations by the "isothermal hold- 
ing" technique, with increasing film thickness the twin 
plates grow in number, and consequently increases the 
hep portion. And again this increase of the hep portion 
occurs not through the conventional martensitic trans- 
formation as in bulk materials, but as the result of an 
increase in the number of twin plates which are shifted 
relative to each other in the [112] direction. Although the 
transition takes place in a rather short time interval of 
10 - 20 ps, it is realized by shear deformation. The final 
structure of the film with a thickness N z = 2Au.c. after 
transformation at T = 430!^ consists almost entirely of a 
hep phase being a system of twin plates. The surface of 
the film obtained has a stepped relief with a step width 
of 2 or 3 atoms. 

In Fig|6]are plotted the atomic radial distribution func- 
tions for such mixed (fec/hep) phases obtained after 
structural transformations and then cooled down to T 
= 100 K, for films of different thickness N z . The dark 
line on bottom plot shows the ARDF for the perfect fee 
lattice, and on top plot the ARDF for the hep structure. 
It is seen from the figure that as the film thickness grows 
from N z — Ylu.c. to N z — 25u.c, the intensity of peaks 
at distances 5.15A and 6.15A corresponding to the hep 
structure gradually increases, and the peak at 6.3 A (fee 
phase) decreases. For the film with N z — 25u.c. the ratio 
between the peak intensities closely corresponds to the 
ARDF of hep lattice. 
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FIG. 6: ARDF for the final fec/hep mixture in thick films of 
different thickness. 



With further increase in the film thickness to N z — 
2bu.c. and more, the sequence and mechanism of trans- 
formations changes again. In this case the bee— >-hcp tran- 
sition is realized by the mechanism of cooperative atomic 
displacements as in bulk zirconium, and the final film 
structure looks quite different. Instead of a distorted 
film with stepped surface relief we observe, as a rule, a 
structure formed by alternating twin hep regions oriented 
in accordance with long-wave vibrations arising at given 
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dimensions of the original crystallite (in FigfJ] the tem- 
peratures of the formation of wavy twin hep structures 
are marked by triangles). So at "slow cooling" there 
exists a very narrow thickness range N z = 17 — 23w.c. 
corresponding to 6.1 - 8.2 nm in which the bec (001) zir- 
conium film transforms to a nanostructured film of fee 
phase. The thickness values for films with a stable twin 
fee phase obtained when simulating structural transfor- 
mations (6.1-8.2 nm) agree with the results of ab initio 
calculations [l7j of the critical size of zirconium particles 
with energetically favorable fee phase (about 6 nm), as 
well as the experimental data (5-10 nm) [l3| . 

Thus the change in the parameter N z , i.e. the thick- 
ness of the bec zirconium film with (001) surface, leads 
to changes in the sequence of structural transitions, the 
mechanism of the bec — > hep transformation, the tem- 
perature of phase transitions, and the final phase. Fur- 
thermore we have found that also the way of setting the 
temperature regime substantially affects the process and 
final structures of phase transformations in films, which 
was not observed in simulation of the phase transitions 
in bulk zirconium with cyclic boundary conditions when 
changing the size of the basic crystallite [29]. 



STABILITY OF FCC ZIRCONIUM FILMS. 
ELASTIC MODULI 

In paper |l| it has been shown that whereas in bulk 
samples the bec lattice is stable, in thin films the fee 
structure may form, since a more closely packed phase 
has a lower surface energy. As stated above, in zirconium 
nanofilms a stable fee phase is realized in a certain tem- 
perature range. Calculations of the total energy in the 
ground state show that fee zirconium has a higher energy 
than the hep phase but lower than the bec lattice. At 
atmospheric pressure in bulk zirconium there exist only 
two phases: hep at low temperatures and bec at high 
temperatures. As shown in Ref.[31|, stabilization of the 
high-temperature bec phase is determined by the con- 
tributions to the free energy from strongly anharmonic 
lattice vibrations with wave vector k = 2/3(111) of the 
Brillouin zone. We do not know any ab initio calculations 
of the ground-state energy in bec and fee zirconium films. 
The total energy of Zr films calculated as a function of 
thickness using the chosen EAM potential at T = K is 
displayed in FigJT] Based on these data one can estimate 
the energy difference A.£/^, cc _^ cc — 

for zirconium films and compare it with that for bulk 
samples. The results obtained show that the presence 
of the surface does not change the sign of AE relative to 
the bulk calculation: for the surface AE^£j cc — 50meV, 
and for the bulk AE%£_ fcc 
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21meV. 

It should be noted that with increasing temperature 
the difference in energy between the bec and fee phases 
decreases : AE^l fcc = 37meV at T = HOOK. (At 



FIG. 7: Total energy of different Zr phases as a function of 
the film thickness at T = K . 



= 1100 the fee phase is metastable, see [29|). Keep- 
ing in mind that the hep zirconium films at T = K are 
energe tically preferable {AE s f u jJ_ hcp = 21meV), one can 
suppose that the fee phase stability is also related to the 
peculiarities of lattice vibrations. So the formation of a 
stable fee phase in Zr films may be connected with the 
contributions to the free energy not only from the sur- 
face states but also from those determined by the lattice 
anharmonicity. 

The stability of fee zirconium is also demonstrated by 
the elastic moduli calculated by various methods. For 
comparison in Table 1 are listed the data for bec and fee 
zirconium structures obtained in our and other papers, 
including the experimental values of the lattice parame- 
ter of fee films. In the second column of the Table are 
indicated the abbreviated names of the calculation meth- 
ods used in the corresponding papers, namely, the ab ini- 
tio calculations WASP, FPLMTO, FPLAPW, CASTEP, 
LDA, and the MD ones using the potentials constructed 
in the embedded-atom model: EAM1 and EAM2. 

Note that in the LDA calculations @ the bec lattice 
was found to be stable at T — K. Our calculations 
show that the bec lattice has a negative shear modulus 
C at zero temperature, i.e. at low temperatures the bec 
lattice is unstable as it must. The experimental data for 
bec zirconium [28| were obtained at T = 915 K. The 
lattice parameter for both the bulk and film samples was 
determined from the total energy minimum in a non- 
relaxed crystallite at T = K. The lattice parameter for 
the fee film was obtained for a non-relaxed film with a 
thickness N z = 10 u.c. 

It should be mentioned here that in films the lattice pa- 
rameter depends on both the film thickness and the tem- 
perature. For example, in the same fee film of thickness 
N z = 10 u.c. as the temperature grows from T = 100 K 
to 500 K, the volume per atom changes from 23.42A 3 to 
23.56A 3 , while in a film with N z — \7u.c. under the same 
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ao 


V 


B 


C 


Cn 


Cl2 


C44 


bcc-bulk 


our calc 


3.563 


22.616 


106.5 


- 9.2 


95.5 


112 


41.9 


WASP [25] 


3.57 


22.75 


89.3 


- 5.5 


82 


93 


29 


LDA [6] 


3.47 


20.87 


93.2 


11.65 


108.7 


85.4 


10.6 


FPLMTO [32] 


3.568 


- 


90.4 


- 


- 


- 


- 


EAM1 [33] 


3.586 


- 


92.4 


- 


- 


- 


- 


EXP [28] 


3.574 


22.83 


96.67 


6 


104 


93 


38 


fcc-bulk 


our calc 


4.544 


23.463 


87.7 


16.53 


109.74 


76.68 


64.33 


FPLAPW [14] 


4.52 


23.1 


91 


21 


119 


77 


53 


EAM1 [33] 


4.542 


- 


87.4 


- 


- 


- 


- 


EAM2 [34] 


4.53 


- 


94 


- 


130 


76 


65 


CASTEP [35] 


4.52 


- 


90.18 


15.495 


110.84 


79.85 


48.24 


fcc-film 


our calc 


4.537 


23.35 


91.2 


15.35 


111.68 


81 


63.47 


EXP [4] 


4.61 














EXP [51 


4.63 














EXP [6] 


4.61 














EXP [71 


5.04 















TABLE I: Table 1. Lattice parameters ao (A), equilibrium volume Vb(A 3 /atom), bulk modulus(B), shear modulus(C') and 
elastic constants (GPa) for Zr in the bcc and fee phases. 





4.68 r 


< 


4.66 


N 


4.64 


O 


4.62 


O 


4.6 - 


in 




S 


4.58 


"S 
E 


4.56 


a 


4.54 


(5 




Q. 


4.52 


CD 




O 


4.5 


i= 


CO 

_l 


4.48 




4.46 



10 



15 20 
N„ unit cells 



25 



30 



FIG. 8: Variation of the lattice parameters (a and c) with the 
fee film thickness at T = 300 K. 



plays the variation of the lattice parameters c (along the z 
direction), and a (along the x, y directions) with increas- 
ing the film thickness. The calculations were performed 
for crystallites relaxed during 100 ps at T — 300 K, and 
then the atomic positions were averaged over the 50 ps 
range for films of different thickness. As seen from the 
figure, with an increase in the film thickness the param- 
eter c diminishes, and the parameter a grows, the ratio 
c/a being equal to 1.042 for the film with N z — 5 u.c, 
and 1.011 for the film with N z — 30 u.c; that is, the 
thinner the film, the greater is the tetragonal distortion 
of the lattice. The calculations of the bulk modulus B, 
and shear constants C and C44 were performed by the 
formulas: 



B = 



d 2 E 



9V da 2 



(1) 



temperature variation the volume changes from 23.45A 3 
to 23.57A 3 . It is seen from the Table that the lattice 
parameters of the experimentally obtained fee zirconium 
films differ from one another, as the films in these works 
were prepared in different ways, under different condi- 
tions, and on different substrates which are known to sub- 
stantially affect the realization and stability of a struc- 
ture, and consequently its lattice parameter. Moreover, 
as shown by our calculations, the bcc and fee films lose 
the cubic symmetry and become tetragonal (bet and fct, 
respectively), the lattice stretches along the z direction 
with free boundary conditions. A tetragonal distortion of 
the film was also found experimentally [6] . Figure [5] dis- 



C" 



d 2 E 



AV da 2 



(2) 



where a is the lattice parameter, ao is the equilibrium 
lattice parameter, V is the volume per atom, E is the 
total energy per atom as a function of the lattice param- 
eter a. The calculation of C44 was also carried out by 
formula (2) but for the (110) orientation of the fee lat- 
tice, i.e. with the x axis along the [110] direction, and 
the y axis along the direction [1-10]. Since the calculation 
was performed for non-relaxed crystallites, the relations 
between the bulk and the elastic moduli correspond to 
the cubic symmetry: 



B 



C*u + 2C12 



C" = 



C n - C 



12 



(3) 



(4) 



In addition, we have calculated the dependence of the 
elastic moduli on the thickness of a non-relaxed fee zirco- 
nium film. Figure |H] presents the size dependences of the 
bulk modulus B, and shear moduli C and C44. The ab- 
scissa is the film thickness N z in unit cells, the straight 
line on the plots corresponds to the bulk value of the 
given quantity. One can see that with increasing film 
thickness the moduli behave in different ways: the bulk 
modulus B diminishes, while the moduli C and C44 in- 
crease, approaching their bulk values. All the moduli 
values are positive, which points to the stability of fee 
films at low temperatures. 



CONCLUSION 

Thus the study performed has evidenced the existence 
of a temperature range (depending on the film thickness) 
in which a stable fee phase may form in nanofilms of pure 
zirconium. The magnitude of tetragonal distortion in this 
phase depends on the film thickness, and ranges from 4% 
to 1% as N z increases from 5 to 30. In films thicker 
than 6.1 nm the fee phase has a twin structure, and a 
hep structure forms at twin boundaries. The portion of 
hep phase increases with the film thickness. At N z > 
32 u.c. (about 12 nm) the fee phase is not observed. An 
analysis of the total energy of different phases in Zr films 
and of the elastic moduli of the bec and fee phases at 
T = K suggests that the fee phase stability at finite 
temperatures is determined by the lattice vibrations. 

Studying the morphology of the structures obtained 
after transformations in bec zirconium films has shown 
that the mechanisms of structural transformations sub- 
stantially depend on the film thickness and orientation. 
For bee (001) zirconium films a limiting thickness of 6.1 
nm was found, at which the mechanism, sequence, and 
temperature of the structural transitions, as well as the 
final phase do change. In thin films the morphology of 
the structure observed after transformations in bec (001) 
films essentially depends also on the film nanostructure 
(in our calculation it was modeled by the basic crystallite 
dimensions in the MD simulation along the x, y axes) . In 
particular, an orthorhombic a" phase was shown to form 
for odd N x>y < 20. 

In thin (up to 6.1 nm) films under slow cooling a "dif- 
fuse" phase transition is observed. The temperature in- 
terval between the onset and the completion of structural 
transition narrows with increasing film thickness, and the 




10 15 20 
N z , u.c. 

FIG. 9: Dependence of the bulk modulus (a), shear modulus 
(b), and C44 (c) on the non-relaxed fee Zr film thickness. 



temperature of the onset of structural transformation de- 
creases. 
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